Counterrotating Dust Disk Around a Schwarzschild Black Hole: 
New Fully Integrated Explicit Exact Solution 
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The first fully integrated explicit exact solution of the Einstein field equations corresponding 
to the superposition of a counterrotating dust disk with a central black hole is presented. The 
obtained solution represents an infinite annular thin disk (a disk with an inner edge) around the 
Schwarszchild black hole. The mass of the disk is finite and the energy-momentum tensor agrees 
with all the energy conditions. Furthermore, the total mass of the disk when the black hole is present 
is less than the total mass of the disk alone. The solution can also be interpreted as describing a 
thin disk made of two counterrotanting dust fluids that are also in agreement with all the energy 
conditions. Additionally, as we will show shortly in a subsequent paper, the above solution is the 
first one of an infinite family of solutions. 
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The observational data supporting the existence of 
black holes at the nucleus of some galaxies, including the 
Milky Way, is today so abundant that there is no doubt 
about the relevance of the study of binary systems com- 
posed by a thin disk surrounding a central black hole. 
Accordingly, a lot of work has been developed in the last 
years in order to obtain a better understanding of the dif- 
ferent aspects involved in the dynamics of these systems 
(see for a recent review of the main works) . Now, due 
to the presence of a black hole, the gravitational fields 
involved are so strong that the proper theoretical frame- 
work to analytically study this configurations is provided 
by the general theory of relativity. Therefore, a strong 
effort has been dedicated to the obtention of exact solu- 
tions of Einstein equations corresponding to thin disklikc 
sources with a central black hole. However, until now, 
any explicit exact solution corresponding to such kind of 
systems has been obtained. 

Stationary and axially symmetric solutions of the Ein- 
stein equations are the best choice to attempt to describe 
the gravitational fields of disks around black holes in an 
exact analytical manner. At the same time, such space- 
times are of obvious astrophysical importance, as they 
describe the exterior of equilibrium configurations of bod- 
ies in rotation. So, through the years, several examples of 
solutions corresponding to black holes or to thin disklike 
sources has been obtained by many different techniques. 
However, due to the nonlinear character of the Einstein 
equations, solutions corresponding to the superposition 
of black holes and thin disks are not so easy to obtain 
and so, until now, exact stationary solutions have not 
been obtained. 

On the other hand, if we only consider static config- 
urations, the line element is characterized only by two 
metric functions. Furthermore, in the vacuum case, the 
Einstein equations system implies that one of the met- 
ric functions satisfies the Laplace equation whereas that 
the other one can be obtained by quadratures. There- 



fore, as a consequence of the linearity of the Laplace 
equation, solutions corresponding to the superposition 
of thin disks and black holes can be, in principle, eas- 
ily obtained. However, only very few solutions has been 
obtained and neither of them has been fully explicitly 
integrated 0, S S H 0, 0, B . 

In this letter we present what seems to be the first fully 
integrated explicit exact solution for the superposition of 
a thin disk and a black hole. We begin by considering 
the Weyl metric for a static axially symmetric spacetime, 
written as lldl 



ds 2 



J*db 2 + e- 2 V dip 2 + e 2A (dr 2 + dz% (1) 



with $ and A only depending on r and z. The Einstein 
vacuum equations leads to the Laplace equation for <f>, 



1 



0. 



(2) 



and, given $, A is obtained by solving the quadrature 

A[$] = J r[($ 2 r - $fjdr + 2§, x $ tZ dz], (3) 

whose integrability is granted by equation (|2j). 
We consider a solution of {2} of the form 

$ = V + <t>, (4) 

where ip corresponds to a black hole solution whereas that 
<\> corresponds to a thin disk solution. Now, by using (|U) 
in ([3]), we obtain 



A[$] = A[V>] + A[4>] + 2A[V>, 



(•5) 



with 



A.[ip,<f>] = J r[(ip jr <t> t r - il>,z<t>,z)dr 

+ (?P,r4>,z + i>, z <f>, x )dz], 



(6) 



2 



a term due to the nonlinear character of ([3]). For 
the black hole we take ip and A[ip] as given by the 
Schwarszchild solution written as 



ib = — In 

V 2 



A[V>1 = -In 
m 2 



c-i 



C + i 
r C 2_! 



(7) 
(8) 



where the prolate spheroidal coordinates are defined by 
means of 



= m 2 (C 2 -l)(l-77 2 ), 



(9) 



and 1 < C < oo, — 1 < rj < 1. 

Now, in order to obtain the thin disk solution, we in- 
troduce the oblate spheroidal coordinates by means of 



r 2 = a 2 (.x 2 + l)(l- y 2 ), 



axy, 



(10) 



with the ranges taken as — oo < x < oo, < y < 1. 
The disk is obtained by taking y = and so is located 
at z = 0, r > a. On crossing the disk, x changes sign 
but does not change in absolute value, so that an even 
function of a; is a continuous function everywhere but has 
a discontinuous x derivative at the disk. By solving the 
Laplace equation (|2|), with the proper boundary condi- 
tions corresponding to a thin disk with an inner edge, we 
obtain for 4> the simple expression 



ay 



a(x 2 + y 2 ) 



(11) 



Then, after a simple integration of ((3]), we obtain for A[cj>] 
the expression 



AM 



a 2 (l -y 2 )A(x,y) 



Aa 2 (x 2 +y 2 ) 4 ' 
where 

A{x, y) = x 4 (9y 2 - 1) + 2x 2 y 2 (y 2 + 3) + y\y 2 



1), 



with a an arbitrary constant and a the inner radius of 
the disk. 

The surface energy-momentum tensor of the disk can 
be written as till 



where V a 
by 



S a .b = eV a V b , (13) 
~*<5". The surface energy density is given 



4a 
a z x° 



Aa 2 x 4 



(14) 



for x > 0, and will be allways positive if we take a > 
0. We then have a dust disk in agreement with all the 




FIG. 1: Energy density e = ae as a function of f = r/a for 
a — a/ a = 1, 9. The first curve on left corresponds to 
5=1, whereas that the last curve on right corresponds to 
a = 9. 



energy conditions. The total mass of the disk can be 
easily computed and we obtain 



M 
2^ 



e(r)rdr = v / 2a^r(l/4), 



(15) 



so that the disk is of infinite extension but with finite 
mass. In Fig. [TJ we plot the energy density for some 
values of a. 

For the combined black hole and disk system, the in- 
tegration of the mixed term A [(b,ip] in §5§ can be done 
with aid of some very useful computational techniques 
introduced in 1121, So, we obtain 



1 



y 



(Ai - A 2 ) , 



(16) 



(12) where 



Ai 



A 2 = 



ay(l - y)(l + x 2 ) - mx(l + y)(( - 1)(1 - rj) 
[ax + m(l — C — v)] 2 + o, 2 (l — y) 2 

ay(l - y)(l + x 2 ) - mx(l + y){( + 1)(1 - rj) 
[ax — m(l + C — v)] 2 + a 2 (l — y) 2 



that vanishes at the axis, when y = 1, and at the disk, 
when y = r) = 0. More details about the obtention of the 
solution will be presented in a subsequent paper. 

The presence of the black hole modifies the energy- 
momentum tensor of the disk, in such a way that arises a 
nonzero pressure at the azimuthal direction. Accordingly, 
the modified energy-momentum tensor of the disk can be 
written as 11 1 



Sab = eV a Vb + pX a X b , 



(17) 



where X a 



The energy density, the azimuthal 
pressure and the "effective Newtonian density" are given, 



3 



respectively, by 



C-i 



C + i 
1 



C + i 

c 



C + i 



(18) 
(19) 
(20) 



where a = e + p. Now, as e > and £ > 1, e and <r 
will be positives everywhere. Accordingly, the energy- 
momentum tensor will be in fully agreement with the 
weak and strong energy conditions. On the other hand, in 
order to fulfill the dominant energy condition, we require 
that p < e, which implies that £ > 2 and a > ^/3m. 
Furthermore, we have that a < e, and thus the total 
mass of the disk when the black hole is present is less 
than the total mass of the disk alone. 

The energy-momentum tensor can also be interpreted 
as the superposition of two counterrotating fluids. In 
order to do this, we cast S ab 



11] 



S ab = e+UlU b + + e-U a _U b _- 



where 



E_ = 



C + i 



r 



(21) 



(22) 



are the energy densities of the two counterrotating fluids. 
The counterrotating velocity vectors are given by 11 1 



Ul = 



ya ± jj X a 

~VT = W' J 



where 



IP 



P 

e 



C-i 



< 1, 



(23) 



(24) 



is the counterrotating tangential velocity. So, we have 
two counterrotating dust fluids with equal energy densi- 
ties. Now, as e± > 0, the two counterrotating dust disks 
are in fully agreement with all the energy conditions. 



As we can see, the above solution presents some very 
interesting properties. First, the associated material 
source presents a very reasonable behavior and its energy 
momentum tensor is in agreement with all the energy 
conditions. Furthermore, its relative simplicity when ex- 
pressed in terms of the spheroidal coordinates, prolates 
and oblates, makes it very easy to study different dy- 
namical aspects, like the motion of particles inside and 
outside the disk, the stability of the orbits and the pos- 
sible existence of singularities. Additionally, as we will 
show shortly in a subsequent paper, the above solution is 
the first one of an infinite family of solutions, all of them 
having many remarkable properties in common with the 
solution here presented. 
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